Modulation spectroscopy with ultracold fermions in an optical lattice 
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We propose an experimental setup of ultracold fermions in an optical lattice to determine the 
pairing gap in a superfluid state and the spin ordering in a Mott-insulating state. The idea is to 
apply a periodic modulation of the lattice potential and to use the thereby induced double occupancy 
to probe the system. We show by full time-dependent calculation using the adaptive time dependent 
density-matrix renormalization group method that the position of the peak in the spectrum of the 
induced double occupancy gives the pairing energy in a superfluid and the interaction energy in a 
Mott-insulator, respectively. In the Mott-insulator we relate the spectral weight of the peak to the 
spin ordering at finite temperature using perturbative calculations. 

PACS numbers: 03.75.Ss 71. 10. Li 71.10.Pm 74.20.Fg 



The use of Feshbach resonances pj and the loading 
of ultracold atoms into optical lattices 0, opened up 
the possibility to realize strong interactions in quantum 
gases. These systems provide outstanding control over 
most of their parameters, which offers the possibility 
to use them as 'quantum simulators' i.e. to inves- 
tigate unresolved questions from quantum many-body 
physics in a well-controllable system. One example is the 
crossover between a Bose-Einstein condensate of bosonic 
molecules and a BCS-like state of fermionic atoms 0,E3- 
A major step in the direction to 'simulate' phenomena in 
periodic potentials has been taken by Kohl et al. [?J and 
very recently by Chin et al. [8j who succeeded in load- 
ing ultracold fermions into a three-dimensional optical 
lattice. 

However, the possibilities to probe ultracold gases in 
optical lattices are still very limited as compared to mea- 
surement techniques in condensed matter systems. On 
top of this some of the known probes for bosonic gases 
are not useful in the context of fermionic gases. Thus 
the lately experimental realization demands for new tech- 
niques. We propose an experimental setup to measure 
the pairing energy of a superfluid state and the inter- 
action strength in a Mott-insulating or liquid state of 
a two component mixture of fermionic atoms in an op- 
tical lattice. Proposals how to gain information about 
the pairing in a superfluid state have been made for 
example in Ref. El [n| for a system without a lat- 
tice and in Ref. |l2j | with a lattice. The setup proposed 
here is based on the experiment by Stoferle et al. 
in which a gas of ultracold bosons was probed using 
the energy absorption induced by periodic modulations 
of the lattice potential. It has been shown that very 
precise information about the excitation structure of a 
bosonic system |H E 113 Eljp, and about the 
degree of incommensurate filling |l9j can be extracted 
by means of this experiment. It would therefore be de- 
sirable to extent the measurement to fermionic systems. 
However, in a fermionic system the energy absorption 
cannot be determined in the same way as it was done 
in the bosonic setup. There the energy absorption im- 
prints a clear characteristic signal in time-of-flight mea- 



surements which is not the case for fermionic systems. 
In the present work we therefore propose to probe a two- 
component fermionic mixture by the observation of the 
double occupancy which is induced by a periodic mod- 
ulation of the lattice potential. The double occupancy 
can be observed by projecting double occupied sites onto 
molecules [2(J ■ We show that this measurement gives in- 
formation on the pair binding energy of the superfluid 
state and the interaction energy in a liquid or Mott- 
insulating state. We further propose how information 
on the spin ordering could be extracted in the case of 
repulsive interaction. 

In our studies we focus on a two component fermionic 
gas which is tightly confined to one dimensional tubes 
along the ^-direction by a strong lattice in the transverse 
directions. Such a system was realized in the experiment 
plf. If an additional periodic potential V x is applied 
along the tubes, the system can be described 0, |23 by 
the Hubbard model 
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Here a — {t, 1} labels the two hyperfine states, cj g and 
Cj.ct are the corresponding creation and annihilation op- 
erators, and hj tCr = ct o .c J a is the number operator on 
site j. The values of hopping parameter J and the in- 
teraction strength U depend among other experimental 
parameters on the lattice potential V x [3l|. The equilib- 
rium properties of the one-dimensional Hubbard model 
have been intensively studied by a variety of techniques 
jUJ Q EE 113 e.g. exact solutions, bosonization, and 
numerical methods. At low temperature in a system 
with attractive interaction a superfluid state is formed, 
whereas for a repulsive interaction the ground state is 
either a half-filled Mott insulator or, away from half- 
filling, a Luttinger liquid. However, the non-equilibrium 
properties are still a challenge. We use the quasi-exact 
adaptive time-dependent density matrix renormalization 
group method (adaptive t-DMRG) H3, IHH^ and per- 
turbative calculations to determine the time-evolution of 
the system. A periodic modulation of the lattice poten- 



2 



tial along the tube V x (t) — Vq (1 + SV cos (cot)), where 
uj is the modulation frequency, can be translated in the 
lattice description to a periodic modulation of J and U, 
i.e. J[T4(£)] and £/[T4(i)]. Thus the system is described 
by an explicitely time-dependent many-body Hamilto- 
nian. In the following Jo and Uq will be used for the 
values J[Vq} and t/[Vb]. 

A system subjected to a periodic modulation of the lat- 
tice potential will on average absorb energy if the modu- 
lation frequency u> corresponds approximately to the en- 
ergy AE needed to create an excitation, i.e. hcu « AE. 
In the limit of strongly attractive interaction one im- 
portant class of excitations is the breaking of pairs. In 
the strongly repulsive limit particle-hole excitations are 
formed. In both cases these excitations cost an energy 
of the order of the interaction energy Uq. Therefore we 
expect that the application of a periodic lattice modula- 
tion can cause energy absorption, if the frequency corre- 
sponds to the excitation energy, i.e. hw s=s Uq. The two 
cases of attractive and repulsive interaction are related 
by an electron-hole transformation for one spin species, 
e.g. c u -» (-rp'cj^., which maps U -> -U HI H|. 
Thus one can restrict the parameters to one sign of the 
interaction only. 
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FIG. 1: Absorbed energy E(t) versus time for different modu- 
lation frequencies for an initial lattice potential height of Vo . If 
the modulation frequency is close to resonance (tku = 1.2E r ), 
energy is absorbed, whereas off resonance (hui = 0.8£" r ) no 
energy absorption takes place. 

In Fig. n the energy absorption for the case of a pe- 
riodic modulation, both close to resonance hui ~ Uq and 
off-resonance is shown. In both cases the absorbed energy 
shows an oscillatory behaviour with the frequency corre- 
sponding to the modulation frequency. Off-resonance on 
average no energy is absorbed by the system, whereas at 
resonance a clear absorption of energy is visible which 
saturates at longer times. 

The energy absorption versus the modulation fre- 
quency is shown in Fig. [21 Here the average value at 
t m = 50h/E r is taken as a measure for the energy ab- 
sorption |32(. For an initial potential Vo = 8E r a clear 
peak centered around the resonance frequency hui w Uq 
is seen (Fig. [2(a)). For a lower initial potential Vq = 4E r 
(Fig. [21 (b)) we see a broadening of the peak and a clear 
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FIG. 2: Energy absorption spectrum (circles) and induced 
double occupancy D (squares) for (a) an initial lattice poten- 
tial Vo = SE r and (b) an initial lattice potential Vo = 4_B r for 
a strong modulation 8V = 0.02 and a s — — 9.2nm. A clear 
peak around hui pa Uo arises in both quantities. For a better 
comparison UqD is shown as well. 

asymmetric form. The maximum of the peak is shifted 
slightly above the value fuo « Uq. These changes are due 
to the fact that the distance of the ground state energy 
to the center of the next band is not exactly given by the 
energy Uq but shows corrections in Jo- 
in the experiment with bosonic atoms [l3l | the broad- 
ening of the characteristic low momentum peak in the 
time-of-flight images has been used as a measure for the 
absorbed energy. In the fermionic experiment the obser- 
vation of the absorbed energy is more difficult. Due to 
Pauli's principle for the momentum distribution all mo- 
menta up to the Fermi momentum are occupied. There- 
fore the absorption of energy will only cause a smear- 
ing out of the step around the Fermi momentum in 
time-of-flight images which can not be used to deter- 
mine quantitatively the amount of energy absorbed |33| . 
We propose to use the double occupancy in the lat- 
tice, i.e. D = ( n j,T n j>i) as a quantity which can 
be observed more easily by projecting pairs of atoms 
onto molecules [20j. Since for strongly attractive (repul- 
sive) interactions the energy absorption corresponds to a 
breaking of pairs (generation of particle-hole excitations) 
we expect that in this limit the energy absorption shows 
up in the expectation value of the double occupancy. In 
perturbative calculations we find that up to the quadratic 
response regime the average induced double occupancy 
and the average absorbed energy are related simply by a 
factor Uq. This result is obtained by expanding J(t) up 
to first order in SV, neglecting the time-dependence in 
the interaction, i.e. U(t) = Uq, and writing the Hamilto- 
nian as H(t) = Hq + g(t)HK, where Hq is the Hubbard 
Hamiltonian at t = 0, Hk is the hopping operator, and 
g(t) = (dJ/dVo)sv=o$V cos(u>t). The linear response to 
this time dependent perturbation gives a purely oscillat- 
ing contribution to the double occupancy, which averages 
to zero. Therefore, we consider the quadratic response, 
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which is given by 



D 2 (t) = — / dt'dt"g(t')g(t") X 2 (t, t', t") (2) 

where \i (*> t'i t") is the quadratic susceptibility. Its rel- 
evant term, i.e. the one that gives a non trivial time de- 
pendence, can be related to the first order susceptibility, 
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Xf{t,t'.t") = -9(t'-t") ± XK (t-t')+ osc. terms, (3) 

where 8(t) is the step function, and Xk{%) = 
(l/ih)([H K (t),H K (0)]). The final expression 

for the double occupancy is = ^2(0) — 

7}{g(0)] 2 ojtImxK(u)/Uo + osc. terms, where Xk{u) 
is the Fourier transform of Xicif)- One finds the same 
linear growth with time for the expectation value of the 
total energy, without the factor 1/Uo. The oscillatory 
terms, however, are different for both quantities. 

In the experimental realization a relatively strong mod- 
ulation strength SV should be used to get an observable 
signal. We checked numerically that even for modulation 
strengths of 20%, the structure of the energy absorption 
spectrum and the spectrum of the induced double occu- 
pancy still contain similar information. In particular we 
find that the position of the peak agrees for both, whereas 
deviations occur in the width and amplitude for higher 
values of the perturbation. Examples for the comparison 
are shown in Fig. [21 

Since we have seen that the experimentally measurable 
quantity, the induced double occupancy, can be used as a 
similar measure as the absorbed energy, we will discuss in 
the following what kind of information can be extracted 
from this quantity. Hereby we focus on (i) how to obtain 
a good measure of the energy necessary to break a pair in 
the superfluid state or to create a particle-hole excitation 
in a Mott-insulator or liquid state and (ii) how to gain 
information about the spin-ordering in a state. 

Direct access to the average energy needed to break 
pairs or to create a particle hole excitations can be ob- 
tained from the position of the absorption peak. At zero 
temperature the dependence of the peak position on the 
initial interaction strength Uo/Jq is shown in Fig. [21 for 
different modulation strengths and fillings (marked by 
squares and cirles) . The results were obtained performing 
the full time-evolution of the system using the adaptive 
t-DMRG [34[ . The position of the peak depends only on 
the initial interaction strength and not on the size of the 
perturbation. For strong interactions the dependence is 
linear. For small values of the interaction strength Uo/Jo 
it deviates from the linear form, since the derealization 
of the fermions becomes more important and the sim- 
ple picture of assuming the state to consist of localized 
particles breaks down |23t |24| . To study the robustness 
of the peak to the influences of finite temperature we 
additionally perform a perturbative calculation for small 



FIG. 3: The position of the peak in the induced double occu- 
pancy. Results obtained by a full time-dependent calculation 
(adaptive t-DMRG) are shown at half filling, a s — 9.2nm for 
a strong perturbation SV = 0.2 (empty squares), at half fill- 
ing for a weak perturbation SV = 0.01 (filled squares), and 
away of half filling for a strong perturbation SV = 0.2 (cir- 
cles). Perturbative results (triangles) are shown for different 
temperatures ksT — 1//3 for small systems sizes (L = 6 and 
N = 6). A good agreement between all results is found, which 
makes the position of the peak a good measure of Uo/Jo- 

system sizes. In Fig. [21 we show the position of the max- 
imum value of the contributions in quadratic response 
calculations [3^ . The collapse of the data for different 
temperatures shows that the position is mainly indepen- 
dent of temperature. The good agreement between our 
results shows that the position is a very robust quantity, 
which is almost independet of temperature, strength of 
the modulation, size and filling of the system. 

In experiments an additional harmonic trapping poten- 
tial is present. We expect that its main influence on the 
induced double occupancy is a broadening of the width 
of the peak, as was found in the bosonic case 0] . In par- 
ticular, the position of the peak will remain unaffected. 
In the limit of strong interaction this can be understood 
by the introduction of potential energy differences be- 
tween neighbouring sites by the trapping potential. This 
causes a shift of the excitation energy for the pair break- 
ing or the creation of particle-hole pairs, but with a sign 
depending on the direction of the excitation. Therefore 
the position of the peak remains unchanged whereas its 
form can be broadened. In the experiment with bosonic 
atoms 13] for strong interaction strength the position 
of the resonance peak in the energy absorption could be 
determined approximately up to an accuracy of 20 Hz. 
Note, that since the broadening of the peak is propor- 
tional to Jo the accuracy is best for strong interaction. 
In the fermionic case it is hard to estimate the loss of ac- 
curacy by the projection onto molecules, but we expect 
that the accuracy which can be obtained is of the same 
order of magnitude. 

In a two component fermionic system with repulsive 
interactions the probability to create a particle-hole ex- 
citation depends strongly on the spin ordering present 
in the intitial state. In a ferromagnetically ordered state 
with only a few domain walls particle-hole excitations are 
strongly supressed by the Pauli principle, whereas in an 
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FIG. 4: The weight of the peak in the induced double occu- 
pancy for different temperatures 1//3 compared to the nearest 
neighbour correlations function iL(rif ,jnu+i) in a small sys- 
tem with L — 6 sites and N = 6 particles. The inset shows 
the competing energy scales influencing the spin ordering. 

antifcrromagnetically ordered state the creation of par- 
ticle hole pairs is facilitated. This means that a state 
which is antiferromagnetically ordered will show more 
absorption than a paramagnetically or ferromagnetically 
ordered state. Thus the weight of the absorption peak 
gives information about the spin ordering. To be more 
specific, for large Uq/Jq perturbative calculations yield 
that the weight is proportional to the nearest neighbour 
correlation (n-ijnij+iS. This can be seen in Fig.0]where 
we compare perturbation results for the weight of the ab- 
sorption and nearest neighbour correlations for different 
temperatures. At large values of Uq/ Jq we find very good 



agreement between the quantities, whereas for smaller in- 
teraction the values deviate. A maximum occurs in both 
quantities at approximately the same finite interaction 
strength. This is due to two competing processes which 
destroy the spin order, the derealization of charge domi- 
nating at low Uo I Jq and the spin exchange dominating at 
large Uq/ Jq. In the inset in Fig.0]this competition of the 
energy scales for the spin ordering, the charge gap and 
the spin exchange coupling are shown. The maximum in 
the weight shows the crossover between the two energy 
scales thus giving information about the spin ordering of 
the system. 

To summarize we have proposed a setup to measure 
the pair-binding energy, the interaction strength, and the 
spin-ordering in a Fermi gas in an optical lattice. The 
proposed experiment observes the double occupancy in- 
duced by periodic lattice modulations via the creation 
of molecules. The robustness of the measurement proce- 
dure to external influences such as temperature or a har- 
monic trapping potential makes an precise experimental 
determination of the pair-binding energy, the interaction 
strength - thereby the scattering length - and the spin 
ordering in the optical lattice possible. Since for strong 
interactions the proposed measurement is based on a lo- 
cal effect, i.e. the creation of particle-hole pairs, we ex- 
pect it to hold as well in higher dimensions. 
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